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ABSTRACT 

Two three-dimensional magnetohydrodynamical simulations of strongly magnetized conical jets, 
one with a poloidal and one with a helical magnetic field, have been performed. In both simulations 
the jet is precessed at the origin to break the symmetry. The jets are denser by an order of magnitude 
than the jets in a previous set of simulations. Theoretical work accompanying the previous simulations 
indicated that an increase in jet density would stabilize the jets beyond the Alfven point and this 
prediction is verified by the present simulations. In the poloidal simulation a significant sheath 
of magnetized moving material developed around the jet and led to additional stabilization. New 
theoretical work analyzing the effect of a jet embedded in a magnetized wind shows that the velocity 
shear, Am = Uj — u e must exceed a "surface" Alfven speed, Vas = [(Pj + pe)(Bj + B^)/(4npjp e )] ll/2 , 
before the jet becomes unstable to helical and higher order modes of jet distortion and this explains 
the enhanced stability found in the poloidal simulation. The fundamental pinch mode is not similarly 
affected and emission knots with spacing ~ 4_Rj developed in the poloidal simulation. Thus, we 
identify a mechanism that can suppress large scale asymmetric structures while allowing axisymmetric 
structures to develop, and astrophysical jets surrounded by outflowing winds will be more stable than 
if surrounded by a stationary or backflowing external medium. Knotty structures along a straight 
protostellar jet like the jet in HH34 or in the inner part of the jet in HH 111 could be triggered by 
pinching of a low magnetosonic Mach number jet surrounded by a suitable wind. Of additional interest 
is the development of magnetic field orientation along the line-of-sight organized by the toroidal flow 
field accompanying helical twisting. On astrophysical jets such structure could lead to a reversal of 
the direction of Faraday rotation in adjacent zones along a jet. Thus, Faraday rotation structure like 
that seen along the 3C465 jet could be attributed to organized magneto-hydrodynamical structures 
produced by the jet flow. 



Subject headings: galaxies: jets — hydrodynamics — instabilities — MHD — ISM: 
jets and outflows 

1. Introduction 

Jet acceleration and collimation schemes imply dynamically strong magnetic fields close to 
the central engine (see Meier, Koide, & Uchida 2001). Numerical studies, e.g., Meier, Payne, & 
Lind (1996), Ouyed, Pudritz, & Stone (1997), Ouyed & Pudritz (1997), Romanova et al. (1997), 
show that the jets created in this fashion pass through slow magnetosonic, Alfvenic and fast 
magnetosonic critical points, whose ultimate velocity may depend on the configuration of the 
magnetic field (Meier et al. 1997) and accelerate up to speeds that may be only a few times the 
Alfven velocity at the Alfven surface (Camenzind 1997). Thus, low magnetosonic Mach numbers 
might be expected for magnetically launched jets (Fendt & Camenzind 1996; Lery & Frank 2000) 
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The numerical studies and also theoretical studies (Begelman & Blandford 1999; Lery et al. 2002) 
indicate that a wind or magnetized wind may surround a more rapidly moving jet. 

All jet flows are susceptible to Kelvin-Helmholtz (KH) or current driven (CD) instabilities. 
The KH instability of 3D jets with purely poloidal or purely toroidal magnetic fields (Ray 1981; 
Ferrari, Trussoni, & Zaninetti 1981; Fiedler & Jones 1984; Bodo et al. 1989) and of jets containing 
force-free helical magnetic fields (Appl & Camenzind 1992; Appl 1996) has been extensively 
investigated for a stationary external medium. Additional investigations have considered the role 
of CD pinching of toroidally magnetized plasma columns (Begelman 1998). The CD instabilities 
are sensitive to the magnetic profile (Appl, Lery, & Baty 2000) and, at least for jets containing 
force-free helical magnetic fields, it appears that the KH instability exhibits faster growth and 
is more likely to be responsible for producing asymmetric structure than CD instability (Appl 
1996) and CD modes appear confined to the jet interior (Appl, Lery, & Baty 2000; Lery, Baty, &; 
Appl 2000). In general, KH growth rates increase as the magnetosonic Mach number decreases 
provided the jet is super-Alfvenic, but the magnetized jet is nearly completely stabilized to the 
KH instability when the jet is sub-Alfvenic (Hardee & Rosen 1999, hereafter Paper I). In Paper 
I the strongly magnetized "light" jets experienced considerable mass entrainment and slowing as 
denser material was entrained following destabilization. Denser jets are predicted to remain stable 
beyond the Alfven point and, in general, have been found to be more robust than their less dense 
counterparts (Hardee, Clarke, & Rosen 1997). 

In this paper we analyze results from two simulations, one with a poloidally magnetized jet 
and the other with a helically magnetized jet. The jets are an order of magnitude more dense 
than those studied in Paper I. In §2 the numerical setup and results of the numerical simulations 
are presented. The simulations are initialized by establishing a cylindrical helically magnetized 
jet across a computational grid in an unmagnetized surrounding medium with pressure gradient 
devised to result in an approximate constant expansion of the jet once pressure equilibrium is 
achieved, i.e., produce a conical jet. In one simulation a significant outwards flowing magnetized 
cocoon develops that influences jet stability and dynamics. This simulation shows how a jet might 
be influenced by a surrounding magnetized wind. Synchrotron intensity images and an integration 
of the line-of-sight magnetic field provide a connection between jet dynamics and observable jet 
emission and polarization structures. In §3 the theory is expanded to include a moving magnetized 
wind medium, and in §4 the predicted effect of density and winds is compared to the simulation 
results. Finally, in §5 we discuss some of the implications for astrophysical jets. 

2. Numerical Simulations 
2.1. Initialization 

Simulations were performed using the three dimensional MHD code ZEUS-3D, an Eulerian 
finite-difference code using the Consistent Method of Characteristics (CMoC) which solves the 
transverse momentum transport and magnetic induction equations simultaneously and in a planar 
split fashion (Clarke 1996). Interpolations were carried out by a second-order accurate monotonic 
upwinded time-centered scheme (van Leer 1977) and a von-Neumann Richtmyer artificial viscosity 
was used to stabilize shocks. The code has been thoroughly tested via MHD test suites as 
described by Stone et al. (1992) and Clarke (1996) to establish the reliability of the techniques. 

All simulations are initialized by establishing a cylindrical jet across a 3D Cartesian grid 
resolved into 130 x 130 x 370 zones. Thirty uniform zones span the initial jet diameter, 2Rq, 
along the transverse Cartesian axes (x-axis and y-axis). Outside the uniform grid zones, the grid 
zones are ratioed where each subsequent zone increases in size by a factor 1.05. Altogether the 
130 zones along the transverse Cartesian axes span a total distance of 30i?o- Along the z-axis 280 
uniform zones span a distance of 40i?o outwards from the jet origin. An additional 90 ratioed 
zones span an additional distance of 40i?o where each subsequent zone increases in size by a factor 
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1.02. The 370 zones along the z-axis span a total distance of 80Rq- Outflow boundary conditions 
are used except where the jet enters the grid where inflow boundary conditions are used. The 
use of a non-uniform grid with larger zones at the grid boundaries has been shown to have the 
beneficial effect of reducing reflections off the grid boundaries as a result of increased dissipation 
of disturbances by the larger zone size at the grid boundaries (Bodo et al. 1995). 

The jets are initialized with a uniform density pj and initial radius Rq. The magnetic field in 
the jet is initialized with a uniform axial component, B z , and a toroidal magnetic component with 
functional form = B^ k sm 2 [irf(r)] where for r < r p k, f(r) = 0.5(r/r p £,) a , and for r p ^ < r < r max , 

f(r) = 1.0 — 0.5 [(1 — r/r max ) /(l — r p k / r max )] b . In these simulations the toroidal component 
increases to a maximum, B vk , at r p ^ = 0.5-Ro, and declines to zero at r max = 0.9i?o so that 
initially all currents flow within the jet. In the simulation with primarily axial field (simulation 
P), a = 1.106 and b = 0.885, and the FWHM of the toroidal field is 0-5.Ro- In the simulation with 
helical magnetic field (simulation T), a = 0.346 and b = 0.277, and the FWHM of the toroidal 
field is 0.8i?o- These particular toroidal profiles are not physically motivated but are well behaved 
numerically and in simulation T the toroidal profile provides a broad cross sectional region within 
the jet in which the toroidal magnetic component is relatively constant (nearly constant helical 
pitch, plasma beta, and magnetosonic speed). In the external medium the magnetic field is equal 
to zero. The equation of hydromagnetic equilibrium 

dr [ P " [ ' + Sir + Sir ) ~ 4„ r ' 

where the term on the right hand side describes the effects of magnetic tension, has been used 
to establish a suitable radial gas pressure profile in the jet by varying the jet temperature. The 
sound, Alfven and magnetosonic speeds are defined as a 2 = Tp/p with T = 5/3, V\ = B 2 /Attp, 
and a^g = a 2 + V\. Since internal dynamics and timescales involve wave propagation across a jet 
with propagation speed a function of jet radius we define radial averages as, for example, 



1 f R 3 

(M ms ) = — / M ms {r)dr , 
tij Jo 



that will differ for different magnetic, density, and temperature profiles. 

In the simulations the external medium is isothermal and the external density p e (z) 
declines to produce a pressure gradient, p e (z) oc p e (z), that leads to a constant expansion, 
Rj(z) = (1 + z/80i?o)i?0) of a constant velocity adiabatic jet containing a uniform poloidal 
magnetic field and an internal toroidal magnetic field that provides some confinement, i.e., 



_ '(Rj/Ro)' 10 / 3 + C p (Rj/R )- 4 - C^Rj/Ro) 

~ [i + c p - cy 



Pe(0) • 



The values of C p and Ca depend on the poloidal and toroidal field strengths, and the ratio of the 
magnetic pressure relative to the thermal pressure. A pseudo-gravitational potential keeps the 
external medium static. Simulation P contains a primarily axial magnetic field configuration with 

only a weak toroidal component, B^ k /B z ~ 0.05 at the inlet, so as to be directly comparable to 
predictions made by a linear stability analysis containing a B z magnetic field component only (§3). 
In simulation T, the axial magnetic field strength is identical to simulation P but now the toroidal 
component B pk /B z ~ 0.44 at the inlet. The jets are initialized with a density pj = 0.20p e (0) and 
with a velocity Uj = 1.78a e that is trans-Alfvenic. The jet sound, Alfven and magnetosonic speeds 
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at jet center normalized to the external sound speed and radial averages of the jet sonic, Alfvenic 
and magnetosonic Mach numbers at the inlet and values of C p and CU for the simulations are 
given in Table 1. Absolute values for jet speed, density, temperature, and magnetic field strength 
are completely determined by choosing values for the external density, p e (0), and the external 
temperature, T e , or sound speed, a e . An absolute length scale is determined by choosing a value 
for i?o- in both simulations the jet is driven by a periodic precession of the jet velocity, Uj, 



Table 1: Initial Conditions 



Simulation 


Uj/a e 


dj/a e 


VA,j/a e 


af s /a e 


(Mj) 


(m/) 




> c p 


Ctf, 


P 


1.78 


1.32 


1.90 


2.31 


1.43 


0.89 


0.75 


1.721 


0.005 


T 


1.78 


1.41 


1.90 


2.37 


1.26 


0.87 


0.70 


1.506 


0.294 



at an angle of 0.01 radian relative to the z-axis with an angular frequency to = 0.5a e /i?o- The 
initial transverse motion imparted to the jet by this precession is well within the linear regime. 
The precession serves to break the symmetry and the precessional frequency is chosen to be below 
the theoretically predicted maximally unstable frequency associated with helical twisting of a KH 
unstable supermagnetosonic jet with jet radius Rj > Rq. In both simulations the precession is in a 
counterclockwise sense when viewed outwards from the inlet. This direction of precession induces 
a helical twist in the same sense as that of the magnetic field helicity. 

In both simulations the jets expand rapidly and after about five dynamical times, 
Td = [a e /Ro)t ~ 5, have achieved an approximate static pressure equilibrium with the surrounding 
medium, and a nearly constant expansion rate on the computational grid with Rj ~ 2Rq when 
z = 80-Rq- The code modifies the initial input static equilibrium state to an appropriate dynamic 




Fig. 1. — Profiles (along the :r-axis) of axial and azimuthal velocity, v z & v y , respectively, density and internal 
energy, p & e, and magnetic field components, B z & B y for simulations P (left) and T (right) at z/Ro = 4.5. The 
vertical lines indicate the location of a conical jet's surface at this distance from the inlet. 
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equilibrium state, i.e., V x (u x B) = throughout the duration of the simulation within an 
axial distance of one jet radius from the inlet. Velocity, density, internal energy and magnetic 
field profiles representative of the initial dynamic equilibrium in the two simulations are shown in 
Figure 1. 

The simulations were terminated at dynamical times (P) = {a e /Ro)t = 96 and (T) 
Td = (a e /Ro)t = 56 when computational difficulties arose. At these times the simulations have 
reached a quasi-steady state out to axial distances of about (P) 60.Ro and (T) 35i?o- At this time 
the jet had undergone (P) ~ 21.4 and (T) ~ 12.5 precessional periods and (P) ~ 3.7 and (T) 
~ 2.2 flow through times, Tfi = (z/u)(a e /Ro), through z = 50i?o- Typically a quasi-steady state is 
achieved after 3 flow through times at a particular distance. Each simulation required roughly 110 
CPU hours on the Cray T90 at the San Diego Supercomputer Center. 

2.2. Simulation Results 

2.2.1. Jet Structure 

Plots of the total pressure, p* = pth + Pbz, normalized by the external pressure p e and velocity 
components along ID cuts parallel to the jet axis in the x-z slice plane are shown in Figure 
2. In this slice plane v x and v y correspond to a radial velocity, v r , and an azimuthal velocity, 
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Fig. 2. — Cuts in the x-z plane parallel to the z-axis at x/Ro = (solid line) 0.233, (dotted line) 0.433, (dashed line) 
0.633, (dot-dash line) 0.833 of jet total pressure, p* , normalized by the external thermal pressure, p e , and the three 
velocity components normalized by the external sound speed, a e , for simulations P (left) and T (right). 
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v<f,, respectively, provided the jet center is not displaced off the z-axis. The jet accelerates in 
response to thermal and magnetic pressure gradients, and to the pseudo-gravitational potential. 
In simulation P with primarily poloidal magnetic field the jet accelerates to 5% above the value 
at the inlet. In simulation T additional acceleration occurs, to about 20% above the inlet value, 
as a result of the "spring" effect from the helical magnetic field. Instability is most obviously 
manifested by the oscillations and growth in the amplitude of the transverse velocity components, 
but also appears in total pressure and axial velocity fluctuations. Low level fluctuations in p* and 
v z appear in both simulations by z ~ 5i?o- Significant fluctuations do not begin to grow until 
z ~ 30i?o in simulation P, but become significant by z ~ 15i?o in simulation T. 

In simulation P we can identify the following features: 

1. At z/Rq < 30 a short wavelength, A <^ 5Rq, low amplitude oscillation in p* (also in v z but 
only one oscillation obvious with the velocity scale used in the figure). 

2. At 30 < z/Rq < 50 a region of rapid growth and complex structure, showing rapid 
oscillations in v y (tu) and a dominant oscillation in v x (v r ) with A <^ 8i?o- 

3. At z > 55i?o dominant out of phase transverse velocity oscillations indicative of helical 
twisting with A ~ 10i?o moving with v w ~ 0.75a e ~ 0.40uj where Uj ~ 1.9o e (v w determined 
from sequential simulation frames). 

In simulation T we can identify the following features: 

1. At z < 40i?o a region with growing amplitude oscillations in p* , v z and v x with A <^ 8Rq. 

2. At z > 40i?o dominant out of phase transverse velocity oscillations indicative of helical 
twisting with A ~ 10i?o moving with v w ~ 0.91a e ~ OAluj where Uj ~ 2.2a e . 



Plots of the axial velocity along with the sonic, Alfvenic and magnetosonic speeds along ID 
cuts parallel to the jet axis in the x-z slice plane inside and outside the jet are shown in Figure 3. 
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Fig. 3. — Cuts in the x — z plane parallel to the z-axis inside the jet (top) and outside the jet (bottom) in simulations 
P (left) and T (right) showing the (solid line) axial flow speed, v z , the (dot-dash line) magnetosonic speed, v ms , the 
(dashed line) Alfven speed, Va, and the (dotted line) sound speed, v s . 
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The cut at x/Rq = 1.69 contacts the conical jet's "surface" at z = 55i?o but adequately represents 
conditions immediately outside the jet at z < 40Rq. The cut inside the jet shows that, in both 
simulations, the jets are super-Alfvenic just outside the inlet and are super-magnetosonic by 
z = 20Rq. The cut outside the jet reveals that a significant boundary layer forms outside the jet 
in simulation P by z = 15Rq with flow and Alfven speeds almost half their values inside the jet. 

The differences in transverse structure in the simulations is shown by velocity, density, 
internal energy and magnetic field profiles at z = 18i?o shown in Figure 4. In both simulations the 
boundary layer has a thickness greater than 0.5i?o but only in simulation P is there significant 
flow and magnetic field in the boundary layer. 



K,6 
K 

1 

>" D.S 

a 

-0.5 

□.s 

□ .4 

m H3 
<\ 0,3 

□ ,1 


0.4 
0,4 
0.S 




1 1 1 1 1 1 1 1 

: p 

; i 


, | , i i , | , , i ij i 

/ ^ 


1 1 1 | 1 1 1 1 


■ i i i 1 i i i 


1 1 1 1 1 1 1 1 1 1 1 1 1 





1 1 1 1 | 1 1 1 

_ \ 

'- _, Tl 


1 1 1 1 1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 1 1 1 1 


■ 1 1 1 II I 1 _ 
/ - 


I 1 

-1 1 1 1 1 1 1 1 





-B 



■ 





1 1 1 1 1 1 1 1 1 1 1 1 1 

B 

/ 


1 1 1 | 1 1 1 l_ 

K ' 


-I I I 1 1 1 1 1 


l 1 l l l l 1 l l l l 1 l 





-3 -E 



Fig. 4. — Profiles (along the a;-axis) of axial and azimuthal velocity, v 




respectively, density and internal 



energy, p & e, and magnetic field components, B z & B y for simulations P (left) and T (right)at z/Ro 
vertical lines indicate the location of a conical jet's surface at this distance from the inlet. 
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2.2.2. Mass Entrainment 

The mass entrained by the jets and the average velocity of jet plus entrained material in the 
two simulations is plotted in Figure 5. In particular, we define the mass per unit length, cr(z), at 
any point along the jet as a = J A fpdydx, where A is the cross sectional area of the computational 
domain at axial position z, and / is set to 1 if the local magnetic field is above 4% of the expected 
maximum field strength along the jet at z [B(z) > 0.04Bj imax (z)] and / is set to otherwise. 
We define the entrained mass by the presence of a magnetic field since only the jet material is 
initially magnetized. The setting of / to 1 or effectively assumes that zones with a fraction of 
jet material, as defined by the presence of magnetic field, are considered mixed with the external 
medium in that zone. While one may expect that flux- freezing will prevent the magnetized mass 
from increasing along the grid, here the jets are unstable. Vorticies at the jet's surface lead to 
intertwining of filaments or sheets of jet and external material, and numerical diffusion leads to 
mixing. This technique provides results similar to estimating mass entrainment by using an axial 
velocity threshold (Rosen, Hardee, Clarke, & Johnson 1999; hereafter RHCJ). Setting the switch 
at a magnetic field strength of 4% of the expected maximum strength at z reduces the sensitivity 
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of the value of a to a small numerical diffusion of the field into the denser "unmixed" external 
medium. 




TO 



Fig. 5. — Profiles of the total magnetized mass (jet plus entrained material) per unit length, a(z), and of the 
average speed of magnetized jet plus entrained material. 



In the very light jet simulations presented in paper I, there was a large increase in a near the 
jet inlet and cr/(Jjet = 10 (&j e t is the expected value at the inlet) was used as a baseline value. 
In these denser jet simulations some immediate increase is observed but now cr/djet ~ 2.5 serves 
as an appropriate baseline value. In simulation P an increase in a/djet and a decrease in {v z )/a e 
at z ~ 30i?o reflects the point at which significant velocity and pressure fluctuations appear. 
Simulation T shows a small increase in a /a jet and a decrease in (v z )/a e at z ~ 10.Ro coincident 
with growth of significant velocity and pressure fluctuations. 

The maximum value of the mass per unit length is (c/crjet) max ~ (P) 4 and (T) 3.5. These 
values imply a jet plus entrained mass of about (P) 1.6 and (T) 1.4 times the initial baseline 
mass per unit length. The mass entrained in simulation P is similar to the poloidally magnetized 
simulation B in paper I. The mass entrained in simulation T is roughly half that in the toroidal 
simulations C & D in paper I, but this difference is attributable to the early termination of 
simulation T before a quasi-steady state could be achieved at large z, e.g., note the decline in 
the jet plus entrained mass at larger axial distances. While the effect of mass entrainment on the 
jets here is somewhat reduced relative to the "light" jets studied in paper I, we again find that 
significant mass entrainment does not occur until the jets destabilize. 

The average axial velocity of magnetized material indicates that the entrained material is 
moving more rapidly than was found in the light jet simulations in paper I. This result indicates 
that a denser jet is capable of accelerating entrained material to higher velocities. In simulation 
P we note that the average axial velocity of magnetized material increases significantly to a 
maximum at z/Rq ~ 25 while at the same time cr/crjet decreases to a minimum. This behavior 
is coincident with the formation of a significant sheath of magnetized material moving at up to 
38% of the jet speed (see Figure 3 for velocities in the sheath). In both simulations the majority 
of the momentum flux, av^, [(P) 98% and (T) 92%] evaluated within the quasi-steady state 
region of the computational grid is carried by jet plus entrained material, i.e., material with 
B(z) > 0.04B jtmax (z). 
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2.2.3. Velocity, Emission & Magnetic Field Images 

Jet axial velocity cross sections shown in Figure 6 illustrate development of the surface 
distortions that promote mixing and mass entrainment, and the large scale distortions that can 
move the jet flow completely off the initial axis. In simulation P small scale surface corrugations 
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Fig. 6. — Jet velocity cross sections at axial distances indicated in the panels for simulation P (top 10 panels) and 
simulation T (bottom ten panels). The number in the upper left corner of each panel indicates distance from the 
inlet in units of in- 
form before the jet becomes grossly unstable. A similar result was obtained in paper I for "light" 
poloidally magnetized jets. An elliptical distortion of the jet cross section is evident in the panels 
at z/Rq = 30 & 35, and considerable cross section distortions and displacement of the jet off the 
initial axis occur at larger z. In simulation T the jet exhibits much reduced distortion and the 
higher order surface corrugations that appear in simulation P at z < 35i?o are suppressed. These 
results are similar to what was found in paper I. An elliptical distortion is evident in the panels 
z/Rq = 10 — 30 and the jet is displaced off the initial axis for z > 25i?o- 

Synchrotron intensity images and images of an integration of the magnetic field component, 
By, along the line-of-sight are shown in Figures 7 &: 8. To some extent these images also reveal the 
extent of jet spreading as only the jet fluid is magnetized. To generate the synchrotron intensity 
images a synchrotron emissivity is defined by Pj(B sin#) 3 / 2 where 6 is the angle made by the 
magnetic field with respect to the line of sight. This emissivity mimics synchrotron emission from 
a system in which the energy and number densities of the relativistic particles are proportional 
to the energy and number densities of the thermal fluid. This simplistic assumption is necessary 
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when the relativistic particles are not explicitly tracked (Clarke, Norman, & Burns 1989). 




Fig. 7. — Synchrotron intensity images in the x-z plane from line-of-sight integrations through the computational 
domain along the t/-axis. Image dimensions are 207?o X 80i?o and the logarithmic scale spans six orders of magnitude. 
Simulations P and T are top and bottom images, respectively. 



Fig. 8. — Integration showing the line-of-sight direction of the magnetic field (B y ) in simulations P (top) and T 
(bottom). Image dimensions are identical to Figure 7. 

The line-of-sight magnetic field changes direction in adjacent bands in simulation P. The 
effect is most pronounced at larger z where a pronounced helical twist has developed in the jet 
flow. The poloidal magnetic field is organized by the toroidal flow field accompanying the helical 
twist. The line-of-sight magnetic image for simulation T reflects the orientation of the initial 
toroidal field component. In simulation T the toroidal magnetic field is too strong to be overcome 
by reorientation of the poloidal field by the toroidal flow field accompanying the helical twist. 
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2.2.4- Pinch Structure at Intermediate Time 

At dynamical times ~ 39 the simulations have not yet developed strong asymmetries. In 
Figure 9 synchrotron intensity images, containing fractional polarization B-vectors formed from 
the Stokes parameters, reveal observational structures that develop on the trans- Alfvenic jets 
at intermediate times. Axisymmetric intensity knots appear contained within the jet interior 




Fig. 9. — Synchrotron intensity images in the x-z plane from line-of-sight integrations through the computational 
domain along the t/-axis. Vectors indicate apparent magnetic field direction. Image dimensions are 10-Ro x 40-Ro and 
simulations P and T are top and bottom images, respectively. 

in simulation P and are closely spaced. Intensity knots do not appear in simulation T and the 
development of sinusoidal oscillation is already apparent in the jet interior. In simulation P the 
emission knots, visible in Figure 9 between 25 < z/Rq < 35, have spacing A ^ 4Ro and move 
with v w ~ 0.85a e ~ 0.45uj where Uj ~ 1.9a e . It is interesting that the emission knots appear in 
simulation P containing only a very weak toroidal field component. In simulation T the stronger 
toroidal field component might have been expected to be more conducive to knot formation via 
a pinch instability. We will show that the presence of flow outside the jet in simulation P has 
suppressed asymmetric structures and allowed the knots to form. 

3. Theoretical Analysis 
3.1. Stability Theory 

Let us model the jet as a cylinder of radius Rj, having a uniform density, pj, a uniform internal 
axial magnetic field, Bj, and a uniform velocity, Uj. The external medium is assumed to have a 
uniform density, p e , a uniform external axial magnetic field, B e , and a uniform velocity, u e . The 
jet is established in static total pressure balance with the external medium where the total static 
pressure is p* = pj + Bj /8tt = p* e = p e + B^/8tt. Inclusion of a small non- uniform toroidal magnetic 
field component like that used in simulation P will not significantly modify results obtained from 
a linear analysis involving small fluctuations incorporating only axial magnetic fields, although 
we expect significant effects associated with the much stronger toroidal magnetic field used in 
simulation T, e.g., Appl &: Camenzind (1992), Appl (1996), Rosen et al. (1999). The general 
approach to analyzing the stability properties of this system is to linearize the one-fluid MHD 
equations of continuity and momentum along with an equation of state where the density, velocity, 
pressure, and magnetic field are written as p = po + p\, u = u + ui p = p + pi, B = B + E$i, 
and subscript 1 refers to a perturbation to the equilibrium quantity. Axial magnetic fields in 



- 12 - 



cylindrical geometry have been investigated in several articles (Ray 1981; Ferrari, Trussoni, & 
Zaninetti 1981; Fiedler & Jones 1984; Bodo et al. 1989) although in general it has been assumed 
that there is no flow in the external medium. The relatively slow jet expansion in the numerical 
simulations is not expected to significantly modify local results based on a completely uniform 
external medium (Hardee 1984). 

In cylindrical geometry a random perturbation of pi, ui 5 p±, and Bi to an initial equilibrium 
state po, uq po, and Bq can be considered to consist of Fourier components of the form 



h(r, 



fi(r) exp[i(kz ±n<p — u>t)] 



where the flow is along the z-axis, and r is in the radial direction with the flow bounded by r = Rj. 
In cylindrical geometry k is the longitudinal wavenumber, n is an integer azimuthal wavenumber, 
for n > the wavefronts are at an angle to the flow direction, the angle of the wavevector 
relative to the flow direction is 9 = tan(n/kRj), and +n and — n refer to wave propagation in the 
clockwise and counterclockwise sense, respectively, when viewed outwards along the flow direction. 
In equation (1) n = 0, 1, 2, 3, 4, etc. correspond to pinching, helical, elliptical, triangular, 
rectangular, etc. normal mode distortions of the jet, respectively. Propagation and growth or 
damping of the Fourier components is described by the dispersion relation 



Pj J±n(Pj R ) _ Pe H±l (PeR) 
Xj J±n(PjR) Xe H { ± ] n (P e R) 



(2) 



where J± n and H±^ are Bessel and Hankel functions, the primes denote derivatives of the Bessel 
and Hankel functions with respect to their arguments, and we have dropped the subscript on the 
jet radius, i.e., R = Rj. In equation (2) 



Pj 



-k 2 + 



{u - ku 3 Y 



{u _ kuj) 2 {a 2 + V 2 Aj) _ k 2 V 2 ja 



1/2 



,Pe 



-k 2 + 



(to — ku e ) A 



(u>-ku e nal + Vl)-k*Vla 



1/2 



and 



Xj = pj[(u - kuj) - k V Aj ] , Xe = p e [(w - ku e y 



k 2 Vl] 



where Uj >e , aj >e and VAj, e are the flow, sound and Alfven speeds in the appropriate medium. 

In general, each normal mode, n, contains a single "surface" wave and multiple "body" wave 
solutions that satisfy the dispersion relation. The behavior of the solutions can be investigated 
analytically and we extend previous analytical 3D results here to include for the possibility of an 
external medium with axial magnetic field and flow. 



3.2. Analytical Approximations 

The n = pinch mode "fundamental" (not a surface wave in this case) wave solution can 
be found from the dispersion relation in the low frequency limit co — > and k — > 0. In this limit 
PjR — > 0, (3 e R — > and the dispersion relation becomes 

Xj « - X e(PjR) 2 [ln(P e R) - m/2]/2 , 

where we have used {j' {l3jR)]/{.h{f3jR)\ -» ~PjR/2 and [H^ (f3 e R)}/[H^ Y (f3 e R)] -> 
(P e R)[ln(P e R) - in/2]. As (kR) 2 -> faster than ln(j3 e R) -oo, the real part of the 
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solution to the dispersion relation remains nearly unmodified by the presence of a magnetized 
wind around the jet and is given by (Paper I) 



V 2 n 2 
V Aj a j 



msj 



± 2 



a 



msj 



msj 



1/2 ^ 



1/2 



(3) 



where a r , 



, ms3 = (a 2 + Vjlj^ ^ is the magnetosonic speed. The imaginary part of the solution is 

vanishingly small in the low frequency limit. These solutions are related to fast (+) and slow 
(— ) magnetosonic waves propagating with (v,j+) and against (uj— ) the jet flow speed Uj, but 
modified by the jet-external medium interface. The unstable growing solution is associated with 
the backwards moving (in the jet fluid reference frame) slow magnetosonic wave. The growth rate 
can only be determined by numerical solution of the dispersion relation. 

The n > higher order mode "surface" wave solutions can also be found from the dispersion 
relation in the low frequency limit. In this case the dispersion relation becomes 



Xj 



-Xe , 



where we have used [j' ±n {P j R)H { ± ) n ((3 e R)]/lJ± n (f3 j R)H i ±n {PeR)} -► -((3 e R)/((3jR). The solution 
is modified by the presence of a magnetized wind around the jet and is given by 



k 



u e + rjUj . 7] 
±i 



1/2 



l + rj 



l + T] 



{Uj - lie) 2 



V 2 
V As 



1/2 



(4) 



where r\ = pj/p e and V\ 8 = (pj + p e )(B 2 + B 2 )/(4irpjp e ). Growth corresponds to the plus 
sign in equation (4), but these higher order surface modes are predicted to be stable when 
(uj — u e ) 2 — V\ s < 0. When there is no flow and no axial magnetic field outside the jet, equation 
(4) reduces to equation (4a) in Paper I. In this case the stability condition can be written as 
uj — (1 +rj)V2j < 0. These surface modes are strongly influenced by flow and magnetic field in the 



medium outside the jet. In the dense jet limit, i.e., rj 



oo and V\ s 



V^Aj^ equation (4) becomes 



u/k f» Uj ± Vaj ■ The present reanalysis in this limit reveals an error in paper I where we stated 
incorrectly that u/k « uj ± Vaj applied when the jet was sub-Alfvenic and stable independent of 
the jet density. The unstable growing solution is associated with the backwards moving (in the 
jet fluid reference frame) Alfven wave. When a jet is trans-magnetosonic and trans-Alfvenic the 
growth (damping) rates of these wave modes can only be determined by numerical solution of the 
dispersion relation. 

The "body" wave solutions are somewhat modified by the presence of a magnetized wind 
surrounding a jet and in the limit u> — > and k ^ the dispersion relation becomes 



where we have used J n (z) p 
a correction term given by 



PjR - nvr/2 - vr/4 



2 ' 



(2/ttz) 1/2 cos(z - nir/2 - 7r/4), 9 = cos~ l e w ir/2 - e, and C n « 1 is 



°™ ~ - - - r(l)' 



* Xj Pe Hk 1 ' (J3 e R) V 2 



TTpjR 



1/2 



f n (f3jR) • 
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The solutions are given by 

B. (" + 2,n + l/2W2- T /2(l-C„) p) 
- (M^/MjMAjY]} - l] 

where m > 1 is an integer When the external medium is unmagnetized and there is no wind 
Xe\uj=o = Pe[u e 2 — Va^ 2 = and C n = 0. Previously we found that unstable body wave 
solutions exist only when the denominator in equation (5) is real and that body mode growth 
rates, with the exception of the 1st pinch body mode, remain small unless the jet is sufficiently 
supermagnetosonic. Thus, we do not expect to find structure attributable to higher order, n > 0, 
body modes in the present simulations. 

3.3. Numerical Solutions 

We have investigated the effects of external magnetic fields and winds on jet stability 
by solving the dispersion relation numerically using root finding techniques. The analytical 
expressions provide an initial estimate for solutions in the low frequency limit. We have considered 
a set of parameters relevant to numerical simulations P & T at a distance z ~ 18i?o from the inlet. 
These parameters are given in Table 2. Solutions to the dispersion relation can be sensitive to 

Table 2: Parameters at z ~ 18i?o 



Simulation 


V 




u e /a™ 


M j 


Mf 


M ms 


Pj/P*i 


Pe/Pt 


P 


0.3 


2.20 


0.77 


2.23 


1.43 


1.20 


0.33 


0.80 


T 


0.3 


2.20 


0.00 


2.23 


1.43 


1.20 


0.33 


1.00 



the wind speed and to the presence or absence of a magnetic field in the wind/external medium. 
Solutions that are growing in one parameter regime can be damped in another parameter regime. 
Thus, we have computed and show both growing and damped solutions to the dispersion relation 
for an unmagnetized and a magnetized wind at wind speeds from zero up to those observed in 
simulation P. The root finding techniques were not always capable of converging to a solution at 
all frequencies or of finding a solution unless a very accurate initial estimate could be supplied. 
Solutions generally consist of a growing and damped pair of roots and it was not always possible 
to find both solutions; e.g., the root finder would not converge to the damped fundamental 
pinch mode solution associated with the outwards moving (in the jet fluid reference frame) 
slow magnetosonic wave. Not all of the possible solutions are shown. For example, the fast 
magnetosonic wave solutions associated with the fundamental pinch mode are not included and 
various additional "wind" associated solutions found in the numerical investigation are also not 
shown. In general, these additional solutions are purely real or damped, and will not affect the 
stability properties of the jet. 

3.3.1. The Fundamental Pinch and Asymmetric Surface Modes 

In the absence of a magnetic field in the external environment the behavior of solutions to 
the dispersion relation is shown for a representative sample of wind speeds in Figure 10. Solutions 
corresponding to the 'backwards' moving waves (in the jet fluid frame) are shown in the leftmost 
panels and solutions corresponding to the 'forwards' moving waves (in the jet fluid frame) are 
shown in the rightmost panels. The wind speeds are indicated in the panels in units of the 
magnetosonic speed in the external medium, in this case equal to the sound speed in the external 
medium. The panels for an unmagnetized external medium with no external wind are indicative 
of the expected stability properties of the jet in simulation T. In simulation T there is little wind 
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or magnetic field external to the jet, i.e., u e = and p e /p* e = 1- The higher order (n > 1) surface 
modes grow much faster than the pinch fundamental mode and would be expected to dominate 
the dynamics. 
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Fig. 10. — Solutions to the dispersion relation for surface pinch (p), helical (h), and elliptical (e) wave modes. 
Dotted lines indicate the real part of the wavenumber. Dashed (growing) and dash-dot (damping) lines indicate the 
imaginary part of the wavenumber. The vertical line indicates the location of the precession frequency used in the 
simulations but normalized to the jet radius and velocity at z = 18i?o • In the panels the wind speed is indicated in 
units of a!™ s . 



As the wind speed is increased the growth (damping) rates of helical and higher order surface 
modes are, at first, greatly decreased as predicted analytically by equation (4) at low frequencies 
(compare the top two rows of panels in Figure 10). For these lower wind speeds growth of the 
fundamental pinch mode is enhanced by the presence of an external wind and can greatly exceed 
that of the helical and higher order asymmetric modes. There is a relatively abrupt change in 
behavior of the solutions when the wind speed is about 50% of the sound speed in the external 
medium. At lower wind speeds the 'backwards' moving waves are growing and the 'forwards' 
moving waves are damped. At higher wind speeds the 'forwards' moving helical and higher order 
surface modes move with speed comparable to the fundamental pinch mode and are growing, and 
the 'backwards' moving waves are damped. At these higher wind speeds growth of the pinch and 
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higher order modes are comparable and relatively large, albeit at high frequencies only. However, 
as the wind speed increases further all growth rates diminish. In general, growth of the higher 
order modes is suppressed as the result of an external wind. 

The solutions to the dispersion relation are very sensitive to a magnetized wind and this 
sensitivity to relatively low speed magnetized winds is shown in Figure 11. As in Figure 10 
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Fig. 11. — Surface mode solutions to the dispersion relation for a magnetized external medium at lower wind 
speeds. Notation is identical to that of Figure 10. 



solutions corresponding to the 'backwards' and 'forwards' moving waves (in the jet fluid frame) 
are shown in the left and right columns, respectively. Here the external magnetic field has a value 
similar to that seen in simulation P. Note that the external magnetic field has a minimal effect on 
the growth rates in the absence of a wind (compare the top left panel in Figure 11 with the top left 
panel in Figure 10). Comparison between the leftmost panels in the top and second row in Figure 
11 shows a reduction in helical and elliptical growth rates at low frequencies predicted analytically 
by equation (4). In general, increase in the wind speed stabilizes the higher order surface modes at 
the lower frequencies while leaving the fundamental pinch mode growth rate unchanged. However, 
when the 'forwards' moving helical and higher order surface modes move with speed comparable 
to the fundamental pinch mode they are growing and the maximum growth rate of the pinch and 
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higher order modes is again comparable (see the bottom panels in Figure 11), albeit growth of 
the higher order modes is restricted to high frequencies. We note here that we do not confirm the 
result for the wave speed of the higher order surface modes in the stable region as was stated in 
paper I. The wave speed is given properly by equation (4) in this paper with the previously stated 
result, v w Uj + Vaj, only appropriate in the dense jet limit. The error in paper I resulted from 
an apparent but not true convergence of the root finding routine. 

There is a relatively abrupt break point in behavior of the solutions when the magnetized 
wind speed is about 35% of the magnetosonic speed in the external medium, and solutions at 
higher magnetized wind speeds are shown in Figure 12. Once again the 'backwards' and 'forwards' 
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Fig. 12. — Surface mode solutions to the dispersion relation for a magnetized external medium at higher wind 
speeds. Notation is identical to that of Figures 10 & 11. Note the different vertical scale in the bottom two panels 
(required to show the greatly reduced growth and damping rates). 



moving solutions are shown in the left and right columns, respectively. In general, higher order 
surface mode growth rates are significantly reduced by these higher speed magnetized winds. In 
these magnetized wind cases our root finding techniques were unable to converge to solutions at 
lower frequencies for the pinch fundamental mode at the highest wind speed and in an intermediate 
frequency range where we assume that the helical surface mode will exhibit a narrow instability 
range similar to that exhibited by the elliptical surface mode. Nevertheless, with the exception of 
the highest wind speeds the fundamental pinch mode is the most unstable and should dominate 
the dynamics. The computations in the bottom two panels indicate that the magnetized wind 
observed in simulation P at the termination of the simulation reduces the growth rate of the pinch 
mode by about an order of magnitude, and restricts growth of the higher order surface modes to a 
very narrow frequency range. Clearly the fast magnetized wind in simulation P is responsible for 
the stability of the jet to far beyond the Alfven point. 
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3.3.2. The Body Modes 

The analytical approximations at low frequency indicate that the real part of the body mode 
solutions will be only minimally modified by the presence of an external magnetic field or wind. 
The extent to which a wind modifies the growth rates at higher frequencies is shown in Figure 
13. Even in the absence of a wind, the low magnetosonic Mach numbers guarantee that body 
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Fig. 13. — Pinch, helical and elliptical first body mode solutions to the dispersion relation for an unmagnetized 
(left column) and a magnetized (right column) external wind medium. Notation is identical to that of Figures 10, 
11, & 12. 



mode growth rates, with the exception of the first pinch body mode, will be small. The first pinch 
body mode has a significant growth rate when the wind speed is low. The addition of a wind and, 
in particular, of a magnetized wind reduces these growth rates even further and at the highest 
magnetized wind speed the body modes are damped. 

4. Interpretation of Simulation Results 

In simulations P and T the flow velocity is sub-Alfvenic for the inlet conditions and 
thus enters the computational grid stable to helical and higher order modes of jet distortion, 
although not stable to the fundamental pinch mode. While in both simulations the jet becomes 
super-Alfvenic within a few jet radii of the inlet, both jets would remain stable beyond this point 
even in the absence of an external magnetized wind as a result of their density. For example, in 
simulation P the jet should remain stable to helical and higher order modes until z ~ 10i?o when 

— (1 + rj)Vjlj > 0, where rj = Pj/p e and V^j = [B? /(iirpj)] 1 / 2 . However, at this distance the 
wind speed has already risen sufficiently, u e ~ 0.2a e , so that (uj — u e ) 2 — (1 + f])V\j < and the 
jet is stable to helical and higher order modes even in the absence of a magnetic field in the wind. 
Simulation P remains stable to helical and higher order modes out to z ~ 30i?o where finally 
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(uj — u e ) 2 - V\ s > 0, where Va s = [(Pj + Pe)(Bj + -B 2 )/(47T/9jp e )] 1 / 2 . Simulation T exhibits no such 
magnetized wind and destabilizes to helical and higher order modes when uj — (1 + Tj)V\ - > at 
z~12i? - 

In general, normal mode solutions come in pairs where the slightly faster moving wave 
(forwards moving in the jet fluid frame) is damped and the more slowly moving wave (backwards 
moving in the jet fluid frame) is growing. As the wind speed increases the spatial growth rate 
of helical and higher order surface modes decreases as the velocity shear decreases and the wave 
speed, measured in the stationary observer frame, increases. However, the fundamental pinch 
mode wave speed (in the jet fluid frame a backwards moving slow magnetosonic wave) and growth 
rate remain relatively unchanged as wind speed increases. At some intermediate wind speed, the 
faster moving (forwards moving in the jet fluid frame) helical and higher order surface mode wave 
speed becomes comparable to the speed of the fundamental pinch mode and now these waves are 
growing whereas the slower moving helical and higher order surface modes are damped. Here the 
fundamental pinch and the helical and higher order surface modes are growing with comparable 
growth rates. The intermediate wind speed value depends on the magnetization of the wind 
and occurs at much lower wind speed for a magnetized wind. As wind speeds rise above this 
intermediate speed all growth rates are reduced, but a larger wind speed is required to achieve 
reduction in the fundamental pinch mode growth rate. 

The different dependence of the fundamental pinch and the higher order surface mode 
growth rates to the presence of an external wind explains the stabilization of simulation P to a 
distance well beyond the Alfven point, and also explains the presence of emission knots in the 
jet in simulation P at intermediate simulation times. In particular, the emission knots, visible in 
Figure 9 at distances 25 < z/Rq < 35, with A <J 4i?o ~ 3.5Rj and with v w ~ 0.85a e ~ 0A5uj 
can be identified with the fundamental pinch mode. If we apply the parameters given in Table 

2 for simulation P we find that the observed knot spacing corresponds theoretically to an 
angular frequency ujRj/uj ~ 1. This frequency corresponds to the maximum growth rate for an 
unmagnetized stationary external medium (see the top left panel in Figure 10) and to the lowest 
frequency at which the growth rate is comparable to the maximum growth rate in the presence 
of a low speed magnetized external wind (see the bottom left panel in Figure 11). At these 
intermediate times and at this distance from the inlet the large magnetized winds seen at the 
end of the simulation have not yet developed. The lower magnetized wind speeds at intermediate 
times, u e ~ 0.15 — 0.35a™ s , suppress the helical and higher order modes leaving the fundamental 
pinch mode to dominate and produce the observed emission knots. The absence of a significant 
wind in simulation T results in no suppression of helical and higher order surface mode instability. 
Thus, in simulation T these normal modes dominate the dynamics at all times and have suppressed 
the development of the fundamental pinch mode and accompanying emission knots. 

In order to make some quantitative comparison with structures observed in the numerical 
simulations at large distances from the inlet we have solved the dispersion relation for parameters 
typical to the jet at z/Rq = 35 & 50. The parameters used for these distances are given in Table 

3 and p e /p* e = 1 and a™ s = a e ; recall that Table 2 contains jet parameters for z/Rq = 18. 

Table 3: Parameters at z ~ 35 & 50Rq 



Sim <£ 


a z/R 


V 


Uj/a e 


aj/a e 


UjRj/Uj 


M j 


Mf 


M ms 


Pj/p* 


P (( 


1 35 


0.35 


2.0 


0.93 


0.36 


2.15 


1.60 


1.25 


0.40 


P (( 


1 50 


0.45 


2.0 


0.85 


0.41 


2.35 


2.00 


1.50 


0.46 


T (< 


i 50 


0.45 


2.2 


1.00 


0.37 


2.20 


1.46 


1.22 


0.35 



In simulation P at z/Rq < 30 the short wavelength, A <J 5Rq, low amplitude oscillation in p* and 
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v z is consistent with the weak fundamental pinch mode instability that remains at the highest 
magnetized wind speeds. The jet destabilizes at about z/Rq = 30 and at 30 < z/Rq < 50 exhibits 
a dominant oscillation in v x (v r ) with A <^ 8Rq ~ 5.3Rj, where Rj ~ 1.5i?o a t z/Rq = 40. The 
observed wavelength is consistent with a growing elliptical distortion at frequency uoRj/uj ~ 0.4 
for the relatively low wind speeds, u e /a™ s <^ 0.2 found at z/Rq ~ 40, and an elliptical distortion 
is evident in the cross sections in panels 30 and 35 in Figure 6. At z > 55Rq dominant out of 
phase transverse velocity oscillations indicative of helical twisting with A ~ 10i?o ~ 6i?j moving 
with v w ~ 0.75a e ~ 0.40uj are consistent with precession induced helical twisting for the observed 
low wind speeds of u e /a™ s ~ 0.2. Thus, the initial precession frequency is preserved through the 
instability point. Movement of the jet center off the z-axis appears at z > 40 in the cross section 
panels in Figure 6. 

In simulation T at z < 40i?o a region of growing amplitude oscillations in p* , v z and v x with 
A <; 8Rq can be readily identified with an elliptical distortion seen in cross sections in Figure 
6 from 10 < z/Rq < 30. We find that the observed elliptical distortion is consistent with the 
elliptical surface mode solutions determined from parameters appropriate to simulation T at 
z/Rq = 18 with A ~ 6.5-R-, when ujRj/iij ~ 0.25. At z > 40i?o dominant out of phase transverse 
velocity oscillations indicative of helical twisting with A ~ 10i?o moving with v w ~ OAluj are 
readily identified with the expected precession driven helical surface mode solution as was found 
for simulation P. The simulations indicate and the theoretical solution to the dispersion relation 
confirms that the helical mode is relatively insensitive to the exact values of the Mach number or 
to the wind magnetization and low wind speed found in the simulations at these large distances 
from the inlet. 

In both simulations the magnitude of the toroidal velocity induced by helical twisting declines 
approximately as exp[— 2(r/Rj — 1)] for r > Rj, and the flow field remains well organized out to 
r > 3Rj even though the magnitude of the toroidal velocity declines rapidly. The direction of the 
toroidal velocity field along the line-of-sight oscillates with a sinusoidal wavelength equal to the 
helical wavelength, i.e., ~ 10i?o- This toroidal velocity field orients the poloidal magnetic field 
in simulation P parallel and anti-parallel to the line-of-sight with the banded structure shown in 
Figure 8. In simulation T a similar toroidal velocity field, albeit with magnitude ~ 2/3 of that 
in simulation P, cannot reorient the poloidal field component to overcome the orientation of the 
initial toroidal field component. 

5. Implications 

We have shown that jets can be stabilized to the KH helical and higher order asymmetric 
normal modes provided the velocity shear, Au = Uj — u e , between the jet and the external medium 
is less than a "surface" Alfven speed, Va s = [{Pj + Pe){Bj + B%) / (4tt pj p e )\ 1 ^ 2 ■ Thus, the presence 
of a magnetic field in the external medium and the presence of an outflowing wind in the external 
medium are stabilizing. At relatively high magnetized wind speeds all normal modes are effectively 
stabilized. At more modest wind speeds or when the jet is sub-Alfvenic the fundamental pinch 
mode remains unstable while the helical and higher order modes can be partially or completely 
stabilized. In this case the fundamental pinch mode dominates the dynamics and emission knots 
can be produced in an initially steady flow, albeit the simulations considered here contain an 
initial precessional perturbation. The initially stable poloidal magnetic field simulation shows that 
the small amplitude precessional perturbation is effectively communicated down the jet to the 
point where the jet becomes unstable. We note here that the wave speed in the stable region is 
not as stated in paper I. The wave speed is given properly by equation (4) in this paper with the 
previously stated result, v w w Uj + Vaji only appropriate in the dense jet limit. The precessional 
perturbation couples to growing helical jet distortions in the super-magnetosonic regime which are 
identical for the poloidal and helical magnetic field simulations. The helical twist leads to a very 
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well organized toroidal velocity field within (see Figure 2) and outside of the jet. The toroidal 
velocity field orients any poloidal magnetic field inside and outside the jet along the direction of 
the toroidal velocity. We note that mixing in simulation P has spread the poloidal field out to 
r <^ 2Rj, but the velocity field remains well ordered out to r ~ 3Rj. Thus, a poloidal magnetic 
field in the medium surrounding a jet could be ordered by a helically induced toroidal velocity 
field out to larger distance than seen in simulation P. 

5.1. Protostellar Jets 

Recent work has suggested that protostellar jets may operate in the low magnetosonic Mach 
number regime as low magnetosonic Mach numbers are expected for magnetically launched jets 
(Fendt & Camenzind 1996; Lery & Frank 2000). However, the low magnetosonic Mach number 
regime is typically KH unstable to helical and higher order normal modes. In such a situation, 
helical and other higher order asymmetric distortions would be expected to dominate the dynamics 
and appearance of these jets. This fact has been demonstrated by somewhat higher Mach number 
numerical simulations (Xu, Hardee, & Stone 2000). While there is evidence for helical and other 
asymmetric distortions in protostellar jets, many of these jets contain emission knots along a 
relatively straight jet, e.g., the knotty HH34 jet which lies within a weak CO outflow (Reipurth et 
al. 1986; Chernin &c Masson 1995; Devine et al. 1997). A more complex example is provided by 
the HH111 jet (Reipurth et al. 1997) which shows emission knots in an initially relatively straight 
jet but with an apparent sinusoidal oscillation and staggered bow shocks at larger distance. Here 
the jet is surrounded by a bubble within a larger molecular outflow (Nagar et al. 1997; Cernicharo 
& Reipurth 1996). The HH34 and HH111 jets are shown in Figure 14. 




Fig. 14. — Images of (top) HH34 and (bottom) HH111 (Devine et al. 1997; Reipurth et al. 1997). 

In addition to HH34 and HH111, observations have shown that other protostellar jets are 
embedded in slower and less collimated molecular outflows (see Richer et al. 2000 and references 
therein). These outflows typically surround but do not necessarily directly interact dynamically 
with the more highly collimated jet, i.e., there is no direct observational evidence for a significant 
wind in the material immediately outside the jet. That there might be a significant wind 
environment around a more highly collimated jet is suggested by studies of jet formation and 
collimation from an accretion disk, e.g., Meier, Koide, & Uchida (2001) and Begelman & Blandford 
(1999) in the extragalactic context and Konigl & Pudritz (2000), Contopoulos & Sauty (2001) 
and Lery et al. (2002) in the protostellar context, or by X-wind models in the protostellar context 
(e.g., Shu et al. 2000). While the properties of the wind close to the jet are important, the manner 
in which the wind is produced and its properties far from the jet are not important. The rapid 
exponential decline in velocity fluctuations outside the jet implies an interaction layer of only 
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a few jet radii in thickness, and it does not matter whether a wind is produced as part of the 
entrainment and mixing process (as in the present simulations) or a wind originates in some other 
fashion (see Lery et al. (2002) for some discussion of disk winds and X-winds, and the relation 
between protostellar jets and molecular outflows). 

Our present results suggest that if the magnetosonic Mach number of a protostellar jet is 
sufficiently low, then a low speed magnetized wind with 10% of the jet speed (or significantly 
higher unmagnetized wind speed) could stabilize the jet to helical and higher order modes of 
asymmetric jet distortion while leaving the fundamental pinch mode to grow and provide a 
trigger for knot spacing not too different from what is seen in the protostellar jets. Note that the 
present results indicate that the wind need extend no more than a few jet radii beyond the faster 
highly collimated jet material. Spatial change in the wind speed such as a lower wind speed at 
larger distance from the source would allow helical and higher order modes to grow, triggered by 
precession at the origin. Thus, it is possible that the interaction between a low magnetosonic Mach 
number protostellar jet and an appropriate outflowing wind could trigger knot formation with 
spacing similar to that observed in the HH34 and HH 111 jets without requiring quasi-episodic jet 
injection on short timescales. 

The present numerical and theoretical results do not include the effects of radiative cooling 
that are important on protostellar jets or the effect of much higher jet density relative to the 
external medium. In general, the inclusion of radiative cooling leads to greater instability (Stone 
et al. 1997; Xu et al. 2000). However, previous 2D (Hardee & Stone 1997) and unpublished 
theoretical results indicate that the presence of significant magnetic field reduces the difference 
between the stability properties of radiatively cooled and adiabatic jets. 

5.2. Extragalactic Jets 

The poloidally magnetized jet in simulation P shows surface corrugations even before it 
destabilizes and after destabilization shows significant cross section distortions in addition to 
elliptical distortion and large scale helical twisting. In simulation T the presence of significant 
toroidal magnetic field component suppresses the filamentation seen on the poloidally magnetized 
jet. While the jet in simulation T exhibits elliptical distortion and helical twisting all higher order 
normal modes are suppressed. This effect was observed in paper I and in supermagnetosonic jet 
simulations containing a significant toroidal magnetic field component (RHCJ). No significant mass 
entrainment occurs in the stable region in the poloidally magnetized simulation P but significant 
mass entrainment accompanies the development of elliptical distortion and helical twisting where 
the jet destabilizes. In simulation T, there is also evidence for mass entrainment beginning 
where the jet destabilizes and elliptical distortion appears. Simulation P shows significantly 
less entrainment and slowing of ordered flow than was observed in the comparable "light" jet 
simulation (simulation B in paper I), and this result confirms that denser jets are more robust. 
Simulation T could not be run long enough to make a similar assessment, but it is clear from 
paper I, the results here, and results from other work (RHCJ) that a significant toroidal magnetic 
field helps to reduce entrainment and helps to maintain jet collimation. 

The line-of-sight appearance of the low magnetosonic Mach number poloidally magnetized 
jets both here and in paper I when they become unstable appears plume like. Simulations at 
somewhat higher magnetosonic Mach numbers with primarily poloidal or weak toroidal magnetic 
field (Rosen, Hardee, Clarke, & Johnson 1999; Rosen & Hardee 2000) also appear plume like. 
This set of poloidal and weak toroidal magnetic field simulations seems most representative of the 
lower power FR I type extragalactic jets whose appearance has been argued (see Bicknell 1994, 
1995) to be the result of significant mass entrainment. Surface turbulence from KH instability has 
been proposed previously as a mechanism for the production of the high observed RM (rotation 
measure) associated with some of the FR I radio sources (Bicknell et al. 1990), although other 
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authors have pointed out quantitative difficulties with this mechanism (Taylor & Perley 1993; Ge 
& Owen 1993). Quantitative difficulties take the form of requiring an excessively large magnetic 
field in a thin Faraday layer or require a Faraday layer thicker than, for example, the dimension of 
the RM structure in order to produce the observed RM with more modest magnetic fields. 

The simulation results presented in RHCJ, Rosen & Hardee (2000) and in paper I revealed 
large scale organization in line-of-sight emission image structure and in the accompanying B-field 
vectors associated with helical twisting and mass entrainment. These organized structures 
extended throughout the mixing region and formed a relatively thick layer around a central jet 
spine. This occurred if polarization B-field vectors were aligned with the jet flow or if polarization 
B-field vectors were transverse to the jet flow if the toroidal magnetic fields were sufficiently 
below equipartition. The present work reveals that the ordering is accomplished by the flow field 
associated with helical twisting, should extend to at least a jet diameter beyond the central jet 
spine, and can produce positive and negative RM bands that alternate and lie across a jet. The 
sum total of all the work suggests that ordering by the flow field outside the jet spine requires 
a primarily poloidal magnetic field in the mixing region and/or in the cocoon environment 
surrounding the jet/tail. Our present simulations do not specifically address the issue of poloidal 
field alignment in the mixing region or in the external medium in extragalactic sources but we 
note that the polarization vectors in the 3C 31 jets suggest a poloidal orientation of the magnetic 
field in a sheath even though the vectors suggest toroidal orientation of the magnetic field in the 
jet spine (Laing 1996). 

Recently Eilek & Owen (2002) have considered the high RM structure associated with radio 
source 3C465 in Abell cluster A2634 shown in Figure 15. In particular the west jet/tail of this 




Fig. 15. — Rotation measure distribution in Abell 2634 (3C 465). The colors are chosen to distinguish positive 
(green & light blue) and negative (red & dark blue) values; 1' is 35.2 kpc (from Eilek & Owen 2002). 
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radio source, whose B-field vectors more or less follow the tail, show positive and negative RM 
bands that alternate and lie across the tail. The emission and RM images can be taken to suggest 
organization resulting from helical twisting of the jet/tail flow. The observations suggest a typical 
scale size for the RM bands of ~ 1" or about 0.7 kpc, and Eilek & Owen find that a similarly 
dimensioned thin ~ 0.7 kpc Faraday layer requires an excessively large, ~ 96 fiG and vastly 
over pressured, Pb/p g ~ 420 magnetic field to produce the observed RM. Thus, they conclude that 
the more likely candidate is a thick Faraday screen associated with ordered much weaker cluster 
magnetic fields. However, in 3C465 the jet/tail diameter, as inferred from the transverse extent of 
the radiating material, at the location of the alternating RM bands is <^ 5" or about 3 kpc. If we 
associate this distance with the thickness of a Faraday screen then the magnetic field required to 
produce the observed RM is reduced to ~ 20/iG and Pb/pg ~ 25. If the jet/tail is not in the plane 
of the sky, the line-of-sight through the sheath will be longer. While the required magnetic field is 
still high, we conclude that RM structure in some radio sources may be associated with large scale 
flow field organization of an ambient and/or jet magnetic field in a thick sheath around a central 
jet spine. 

P. Hardee and A. Rosen acknowledge support from the National Science Foundation through 
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